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>n . Abstract 

CO I The perturbative QCD predictions for the small x behaviour of the nucleon 

structure functions F2^l{x,Q'^) and gi{x,Q'^) are summarized. The importance of 
If^ , the double logarithmic terms for the small x behaviour of the spin structure function 

I gi{x, Q'^) is emphasised. These terms correspond to the contributions containing the 

' leading powers of as^n(l/x)^ at each order of the perturbative expansion. In the 

non-singlet case they can be approximately accounted for by the ladder diagrams 
Qh! with quark (antiquark) exchange. We solve the corresponding integral equation 

with the running coupling effects taken into account and present estimate of the 
effective slope controlling the small x behaviour of the non-singlet spin structure 
function gi{x, Q^) of a nucleon. 
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1. Introduction 



Analysis of the deep inelastic scattering structure functions in the limit of small values 
of the Bjorken parameter x allows to test perturbative QCD in a new and hitherto unex- 
plored regime |l|, 0, ||. The advent of the HERA ep collider has opened up an experimental 
possibility to confront the theoretical QCD expectations with the experimental data on 
deep inelastic scattering in the small x region. The relevant theoretical framework is 
provided in this case by the Balitzkij, Fadin, Kuraev, Lipatov (BFKL) equation for the 
unintegrated gluon distribution f{x,Qf), where Qt denotes the transverse momentum of 
the gluon and x its longitudinal momentum fraction [Q. This equation sums the leading 
powers of ln{l/x) at each order of the perturbative expansion i.e. it corresponds to the 
leading ln{l/x) approximation. The small x behaviour of the structure functions is driven 
by the gluon through the g ^ qq transitions ||^. 

Perturbative QCD predicts a strong increase of structure functions with the decreasing 
parameter x and the experimental data from HERA are consistent with this prediction 
0, 0. This increase is much stronger than that which would follow from the expecta- 
tions based on the "soft" pomeron exchange mechanism with the soft pomeron intercept 
ctsoft ~ 1-08 as determined from the phenomenological analysis of total hadronic and real 
photoproduction cross-sections 

The small x behaviour of structure functions for fixed reflects the high energy 
behaviour of the virtual Compton scattering total cross-section with increasing total CM 
energy squared since W'^ = Q'^{l/x — l). The Regge pole exchange picture is there- 
fore quite appropriate for the theoretical description of this behaviour. The high energy 
behaviour which follows from perturbative QCD is often refered to as being related to the 
"hard" pomeron in contrast to the soft pomeron describing the high energy behaviour of 
hadronic and photoproduction cross-sections. 



The pomeron does not contribute to the spin dependent structure function gi{x,Q'^) 
which is controlled by the exchange of reggeons corresponding to axial vector mesons 



I0| . The pomeron also decouples, of course, from the "non-singlet" part of the unpolar- 

2-^2 



ized structure functions i.e. from the combination Ft' — which at small x is controlled 



by the A2 reggeon exchange. The small x behaviour of the structure function gi{x,Q'^) 
in perturbative QCD is sensitive to a new class of double logarithmic ln{l/x) contribu- 
tions i.e. to the terms which contain powers of asl'n?{l/x) in the perturbative expansion 



11, 12, 13, 14, 15 1 . These terms invalidate naive Regge pole model expectations and 



generate the singular small x behaviour of the structure function gi{x, Q 



The purpose of this paper is to briefly summarize some of the QCD expectations con- 
cerning the small x behaviour of the unpolarized and polarized structure functions of the 
nucleon. In the next section we recall the Regge pole model expectations for the small x 
behaviour of both the polarized and unpolarized structure functions while section 3 is de- 
voted to the discussion of the BFKL equation 0] and its generalization based on angular 
ordering, i.e. the Catani, Ciafaloni, Fiorani and Marchesini (CCFM) equation |I6| , |1^ |18[. 
In Sec. 4 we discuss the small x behaviour of the polarized structure function gi^x.Q'^) 
concentrating for simplicity on its non-singlet part. We recall the (approximate) inte- 
gral equation which resums the double logarithmic terms and discuss its solution after 
taking into account the running coupling effects. We present an estimate of the effec- 
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tive slope A controlling the small x behaviour of the spin-dependent structure function 
gi^{x, Q^) = gi{x, Q^) — gi{x, Q^) ~ x^^. Sec. 5 will contain our main conclusions. 

2. Regge pole model expectations for the small x be- 
haviour of structure functions 

The Regge pole model describes the high energy behaviour of scattering amplitudes in 
terms of the exchange of reggeons which correspond to poles in the complex angular mo- 
mentum plane in the crossed channel [^. The total cross-sections are related through 
the optical theorem to the imaginary part of the scattering amplitudes in forward di- 
rection (i.e. for t = where t denotes the momentum transfer squared). The Regge 
pole model description of the scattering amplitudes when combined with the optical the- 
orem gives therefore automatically predictions for the high energy behaviour of the total 
cross-sections: 

I / e \ ai(0)-l 

cr*f (s) ~ -ImAa,,ab{s, ^ = 0) = E A - • (1) 

S ^ \SoJ 

The quantities ai(0) are the intercepts of the Regge trajectories and Pi denote the cou- 
plings. 

The Regge pole exchange is to a large extent a generalization of the single particle 
exchange and so the reggeons carry the same quantum numbers as hadrons which partic- 
ular Regge pole trajectories correspond to. There is an obvious hierarchy of Regge pole 
contributions depending upon the magnitude of their intercept. The contribution with 
highest intercept corresponds to pomeron which carries vacuum quantum numbers. 

The high energy behaviour of the total hadronic and (real) photoproduction cross- 
sections can be economically described by two contributions: an (effective) pomeron with 
its intercept slightly above unity (~ 1.08) and the leading meson Regge trajectories with 
intercept aR^O) ~ 0.5 The reggeons can be identified as corresponding to p,u, f 
or A2 exchange (s) depending upon the quantum numbers involved. All these reggeons 
have approximately the same intercept. One refers to the pomeron obtained from the 
phenomenological analysis of hadronic total cross sections as the "soft" pomeron since 
the bulk of the processes building-up the cross sections are low pt (soft) processes. 

The Regge pole model gives the following parametrization of the deep inelastic scat- 
tering structure function F2{x,Q'^) at small x 

F,{x,Q^) = Y.UQ^)x^-''^^'\ (2) 

i 

The relevant reggeons are those which can couple to two (virtual) photons. The (singlet) 
part of the structure function F2 is controlled at small x by pomeron exchange, while the 
non-singlet part ^2^*^ = F2 ~ ^2 by ^2 reggeon. Neither pomeron nor A2 reggeons couple 
to the spin structure function gi{x, Q"^) which is described at small x by the exchange of 
reggeons corresponding to axial vector mesons i.e. to Ai exchange for the non-singlet 
part gi^ = gi — gi etc. 

grM')=im^-'''^^'^ (3) 
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The reggeons which correspond to axial vector mesons are expected to have very low 
intercept (i.e. < etc.) 

Several of the Regge pole model expectations for structure functions are modified by 
perturbative QCD effects as will be briefly described in the forthcoming Sections. 



3. BFKL pomeron and QCD predictions for the small 
X behaviour of the unpolarized structure function 

At small values of the parameter x the dominant role is played by the gluons. The basic 
quantity is an unintegrated gluon distribution /(x, Q^) where x denotes the momentum 
fraction of a parent hadron carried by a gluon and Qt its transverse momentum. The 
unintegrated distribution f{x,Qf) is related in the following way to the more familiar 
scale dependent gluon distribution g{x,Q'^): 



xg{x,Q')=l ^/(x,Q?). (4) 



Q't 

In the leading logarithmic approximation the unintegrated distribution f{x,Qf) satisfies 
the BFKL equation Q which has the following form: 



as 



1 dx' r <Pq 
X x' J nq^ 



f{x,Ql) = f{x,Ql)+ 

/(x', {q + Q,f) - fix', QlmQl - 



[q + Qt) 



(5) 



where 

Us = . (6) 

71 

The first and the second terms in the right hand side of eq. (j^) correspond to the real 
gluon emission with q being the transverse momentum of the emitted gluon, and to the 
virtual corrections respectively. f^{x,Ql) is a suitably defined inhomogeneous term. 

After resumming the virtual corrections and " unresolvable" gluon emissions (g^ < 
11^) where is the resolution defining the "resolvable" radiation, equation (|^) can be 
rearranged into the following "folded" form: 



/(a:.«?) = /°(a:.«?) + 



a., 

where 

\ Jz Z' JfjP- 

and 
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Equation (|^ sums the ladder diagrams with reggeized gluon exchange along the chain 
with the gluon trajectory adQt) = 1 ~ ^^'^{Q't 1 1^'^) ■ 

For the fixed coupling case eq.(|^) can be solved analytically and the leading behaviour 
of its solution at small x is given by the following expression: 

/(.,0?)~(«),-^e.p(-^J™j (10) 

with 

\bfkl = 4/n(2)a, (11) 

A" = a,28C(3) (12) 

where the Riemann zeta function C(3) ^ 1.202. The parameter Q is of nonperturbative 
origin. 

The quantity 1 + Xbfkl is equal to the intercept of the so called BFKL pomeron. Its 
potentially large magnitude (~ 1.5) should be contrasted with the intercept asoft ~ 1-08 
of the (effective) "soft" pomeron which has been determined from the phenomenological 
analysis of the high energy behaviour of hadronic and photoproduction total cross-sections 



In practice one introduces the running coupling (XsiQl) the BFKL equation (|). 
This requires introduction of the infrared cut-off that would prevent entering the infrared 
region where the coupling becomes large. The effective intercept Xbfkl found by numer- 
ically solving the equation depends on the magnitude of this cut-off [|T^ . 

The structure functions F^^hi.x^Q'^) are described at small x by the diagram of Fig.l 
which gives the following relation between the structure functions and the unintegrated 
distribution /: 

i^2,.(x, Q') =f^4l ^F^:l{x', Ql Q')n^, Ql). (13) 

Jx X J L/f X 



The functions F2f;{x', Qt,Q^) may be regarded as the structure functions of the off-shell 
gluons with virtuality Q^. They are described by the quark box (and crossed box) dia- 
gram contributions to the photon-gluon interaction in the upper part of the diagram of 
Fig. 1. The small x behaviour of the structure functions reflects the small z {z = x/x') 
behaviour of the gluon distribution f{z, Ql). 

The equation ( [TBD is an example of the " kt factorization theorem" which relates mea- 
surable quantities (like DIS structure functions) to the convolution in both longitudinal as 
well as in transverse momenta of the universal gluon distribution /(z, Qf) with the cross- 
section (or structure function) describing the interaction of the "off-shell" gluon with the 
hard probe [^. The kt factorization theorem is a basic tool for calculating observable 
quantities in the small x region in terms of the (unintegrated) gluon distribution / which 
is the solution of the BFKL equation. 

A more general treatment of the gluon ladder is provided by the CCFM equation 



based on angular ordering along the gluon chain IT^, [l^. This equation embodies both 
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the BFKL equation at small x and the conventional Altarelli-Parisi evolution at large x. 
The unintegrated gluon distribution / now acquires dependence upon an additional scale 
Q which specifies the maximal angle of gluon emission. The CCFM equation has a form 
analogous to that of the "folded" BFKL equation (^: 

f{x,QlQ') = f\x,QlQ')+ 

"^r^ / 00W-^-/-OAK(^,Q?,g^)^^^/(x',(g + Q.f,g^)) (14) 

where the theta function Q{Q — qx/x') refiects the angular ordering constraint on the 
emitted gluon. The "non-Sudakov" form-factor A/j(z, , q^) is now given by the following 
formula: 

^R{z,Ql,q^) = exp 

Eq. (|T^ still contains only the singular term of the g ^ gg splitting function at small z 
yet its generalization which would include the remaining parts of this vertex (as well as 
quarks) is possible. 



J z Z J Q 



(15) 



In Fig. 2 we show the results for the structure function F2 calculated from the kt 
factorization theorem with the function / obtained from the CCFM equation |pi[. We 



confront these predictions with the most recent data from the HI and ZEUS collaborations 
at HERA p, 1^ as well as with the results of the analysis which was based on the Altarelli- 
Parisi equation alone without the small x resummation effects being included in the 



formalism ||22| , P3| . In the latter case the singular small x behaviour of the gluon and sea 
quark distributions has to be introduced in a parametrization of the starting distributions 
at the moderately large reference scale = Ql (i.e. Ql ~ AGeV^ or so) [^. One can also 
generate the singular behaviour dynamically starting from the non-singular "valence-like" 
parton distributions at some very low scale /Xq = 0.35GeV^ pSf. In the latter case the 



gluon and sea quark distributions exhibit the following "double logarithmic behaviour" 



xg{x, Q2) ^ exp aQ')ln{^)j (16) 
where the evolution length ^(Q^) is defined as below: 

For very small values of the scale /Xq the evolution length ^(Q^) can become large for 
moderate and large values of and the "double logarithmic" behaviour (0) is within 
the limited region of x, similar to that corresponding to the power like increase of the 
type A ~ 0.3. This explains similarity between the theoretical curves presented in 
Fig.2. 
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4. Small X behaviour of the nonsignlet unpolarized 
and polarized structure functions 

The discussion presented in the previous Section concerned the small x behaviour of the 
singlet structure function which was driven by the gluon through the g ^ qq transition. 
The increase of the gluon distribution in the small x limit implies similar increase of the 
structure function F2{x,Q^). It turns out to be stronger than the increase that would 
follow from the "soft" pomeron exchange with its relatively low intercept asoft ~ 1.08. 

The gluons of course decouple from the non-singlet channel and the mechanism of 
generating the small x behaviour in this case is different. 

The simple Regge pole exchange model predicts in this case that 

F2^^(x,Q2) = Fi{x,Q') - F-ix,Q') ~ (18) 

where ayi2(0) is the intercept of the A2 Regge trajectory. For 0^2(0) ~ 1/2 this behaviour 
is stable against the leading order QCD evolution. This follows from the fact that the 
leading singularity of the moment ■Jqqiuj) of the splitting function Pqq{z): 

7(^) = /' -z-P,,{z) (19) 

JO z 

is located at u = and so the (nonperturbative) A2 Regge pole at u = ^^2(0) ~ 1/2 
remains the leading singularity controlling the small x behaviour of the non-singlet struc- 
ture function. 



The novel feature of the non-singlet channel is the appearence of the double log- 
arithmic terms i.e. powers of as/n^(l/x) at each order of the perturbative expansion. 
These double logarithmic terms are generated by the ladder diagrams with quark (anti- 
quark) exchange along the chain. The ladder diagrams can acquire corrections from the 
" bremsstrahlung" contributions [|l^, 15 1 which do not vanish for the polarized structure 
function g^^x,Q^) 0. 

In the approximation when the leading double logarithmic terms are generated by 
ladder diagrams illustrated in Fig. 3 the unintegrated non-singlet quark distribution 
f^^{x, kf) (g^"^ = u + u — d — d) satisfies the following integral equation : 



where 
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as = —as (21) 
Sn 

and Ql is the infrared cut-off parameter. The unintegrated distribution f^^{x, Qf) is, as 

usual, related in the following way to the scale dependent (nonsinglet) quark distribution 
g^^(a;,Q2): 



dQ] 



g^^(x,g^)=/ ;^/r(x,Q?). (22) 
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The upper limit Ql/z in the integral equation (EDI) follows from the requirement that 



the virtuality of the quark at the end of the chain is dominated hy Qf. A possible non- 
perturbative A2 reggeon contribution has to be introduced in the driving tem i.e. 



(23) 



at small x. 



Equation (pO|) implies the following equation for the moment function f^^{uj,Q 



Ql Q 



r2\ 



(24) 



(25) 



Equation (^4D follows from (^) after taking into account the following relation: 

i' ? (i -')-^ - * " (§)" - ■ 

For fixed coupling equation ( p^ can be solved analytically. Assuming for simplicity 
that the inhomogeneous term is independent of Qf (i.e. that f^^{uj, Ql) = C{lj) ) we get 
the following solution of eq. 



fJ'iu;,Q-t) = Ciuj)Ri&,,uj) 



2\7 ias,t^) 



where 



and 



7 («.,'^) 



00 — \/uj'^ — 4q;s 



u;7 (as,u;J 
a. 



(26) 



(27) 



(2^ 



Equation (|27|) defines the anomalous dimension of the moment of the non-singlet quark 
distribution in which the double logarithmic IniXjx) terms i.e. the powers of ^ have 
been resummed to all orders. It can be seen from ( ^Tf ) that this anomalous dimension has 
a (square root) branch point singularity at = 



(29) 



This singularity will of course be also present in the moment function (lu, Q^) itself. It 
should be noted that in contrast to the BFKL singularity whose position above unity was 
proportional to as, Cu is proportional to ^/a^ - this being the straightforward consequence 
of the fact that equation ( P^ sums double logarithmic terms (%)". This singularity gives 
the following contribution to the non-singlet quark distribution f^^{x,Qf) at small x: 



X 



/n3/2(l/x) 



(30) 



For small values of the QCD coupling this contribution remains non - leading in compar- 
ison to the contribution of the A2 Regge pole. 

As has been mentioned above the corresponding integral equation which resums the 
double logarithmic terms in the spin dependent quark distributions is more complicated 



8 



than the simple ladder equation (^) due to non-vanishing contributions coming from 
bremsstrahlung diagrams. It may however be shown that , at least as far as the non- 
singlet structure function is concerned, these contributions give only relatively small 
correction to uj. In what follows we shall therefore limit ourselves to the simple lad- 
der equation ( pUD assuming that it will describe the spin dependent parton densities 
as well. The inhomogeneous term will now however be different and will contain the 
Ai reggeon contribution. We will limit ourselves to the non-singlet structure function 
g^^{x, Q2) = gli^x, Q2) _ gn{^x^ Q^) = where 

Aq^^ix, Q^) = Au{x, + Au{x, Q2) - Ad{x, Q^) - Ad{x, Q^) (31) 

where Am, Ad and Am, Am denote the corresponding spin dependent quark (antiquark) 
distributions. We will present an estimate of the effective slope of the non-singlet distribu- 
tions after numerically solving equation (pOD taking into account the asymptotic freedom 
corrections i.e. allowing the coupling constant Us to run. 

The main interest in applying the QCD evolution equations to study the spin structure 
function is that the naive Regge pole expectations based on the exchange of low-lying 
Regge trajectories become unstable against the QCD perturbative "corrections". The 
relevant reggeon which contributes to g\^'^{x,Q^) is the Ai exchange which is expected 
to have a very low intercept aAi(O) < 0. The perturbative singularity generated by 
the double logarithmic ln{l/x) resummation can therefore become much more important 
than in the case of the unpolarized case when it is hidden behind leading A2 exchange 



contribution. Even if we restrict ourselves to the leading order QCD evolution [^, |25 
then the non-singular x~°'^i^^^ behaviour (with a^^(O) < ) becomes unstable as well and 
the polarized quark densities acquire singular behaviour: 



Ag^^(x, Q2) ^ exp{2^^^S{Q^)ln{l/x)) (32) 

where 

eHQ') = f^<^s{q') (33) 

This follows from the fact that APqq{z) = Pqq{z) where Pqq{z) and APqq{z) are the split- 
ting functions describing the evolution of spin independent and spin dependent quark 
distributions respectively and from the fact that Pqq{z) — > const as 2; — > 0. 



The introduction of the running coupling effects in equation ( P^ turns the branch 
point singularity into the series of poles which accumulate at a; = 0. If one makes the 
substitution dig as{Qt) then the corresponding equation for the moment can still be 
solved analytically by reducing it to Kummer's differential equation with the bound- 
ary conditions fixed by the requirement that the solution should match the perturbative 



expansion of the original integral equation [12|. It has however been argued that the the- 



oretically more justified introduction of the running coupling is through the substitution 
as — > Q-siQll z) under the integrand on the rhs. of eq. (|20| ) ||14|| . We perform a numerical 
analysis with both prescriptions and estimate the effective slope controlling the small x 
behaviour of the solution. The two prescriptions of introducing the running coupling ef- 
fects into the double ln{l/x) resummation lead to the following integral equation(s) for 
A/f^(a:,Q?) 

A/f -(X, Ql) = A/5^(x, g?) + &M) f - It ^A/r(-' Qt) (34) 
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or 

A/f ^(x, g?) = A/5-(x, g?) + /' ^a.(g?A) ^A/f -(^, gf ). (35) 

We solve these equations numerically, suitably adapting the method developed in P(]| . 
We assume for simplicity that the inhomogeneous term Af^Q^{x,Ql) is independent of 
both X and of Q^. The flat x dependence of the inhomogeneous term corresponds to 
the assumption that ayij(O) = 0. We assume the LO parametrization of the running 
coupling with four flavours and set A = 0.2GeV. The cut-off parameter Ql was set equal 
to Ql = lGeV\ In Fig.4 we show the effective slope \{x, Ql) for Ql = lOGeV^, 

A x,g = ' ) 36 

aln{l/xj 



for the solutions of equation p^ ) and of equation (pSl). We also show the effective slope 
corresponding to the approximation of retaining only the first integral on the right hand 
side of eq. (P^. In this approximation one sums the single logarithmic ln(l/x) terms 
accompanied by the leading powers of C,^^{Qt)- The small x behaviour is then asymp- 
totically given by eq. . We have also solved eq.(P3D assuming that the driving term is 
given by the A2 reggeon contribution (see eq.(|23|)). The resulting slope of the solution is 
also shown in Fig. 4. It can be seen that the second prescription (i.e as —>■ as{Ql/z)) 
for introducing the running coupling effects makes the effective slope smaller than in the 
case when one makes the substitution as — > asiQl). The perturbative Regge singulari- 
ties are not very important for the unpolarized structure functions and the (input) x"^/^ 
behaviour is not altered substantially by the perturbative Regge singularity. The double 
logarithmic resummation is however very important for generating the singular small x 
behaviour of the polarized structure functions. The results of our estimate suggest that a 
reasonable extrapolation of the (non-singlet) polarized quark densities would be to assume 
an x~^ behaviour with A ~ 0.2 0. Similar (or more singular) extrapolations of the spin- 
dependent quark distributions towards the small x region have been assumed in several 



recent parametrizations of parton densities |26, 27, 28, 2£ . 



5. Summary and conclusions 

In this paper we have briefly summarized the theoretical QCD expectations concerning 
the small x behaviour of the deep inelastic scattering structure functions in both unpo- 
larized and polarized deep inelastic scattering. In the latter case we have for simplicity 
focussed on the non-singlet structure functions which, at small x , can be (approximately) 
described by the ladder diagrams with the quark (antiquark) exchange. We have solved 
the corresponding integral equation taking into account asymptotic freedom effects and 
estimated the effective slope controlling the small x behaviour of the non-singlet structure 
function. The perturbative QCD effects become significantly amplified for the singlet spin 
structure function due to the mixing with the gluons. The simple ladder equation may 
not however be applicable for an accurate description of the double logarithmic terms in 
the polarized gluon distribution AG. 



^ The effective slope A turns out to depend weakly upon the infrared cut-off parameter Qg. If we set 
Qo = then A — > (33-^2jv 1 ~ 0.32. I thank Misha Ryskin for pointing out this fact to me. 
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Figure captions 

1. Diagrammatic representation of the kt factorization formula (0). 

2. A comparison of the HERA measurements of -F2 [§|, 0] with the predictions based 
on the kt factorization formula (|13D using for the unintegrated gluon distributions / 
the solutions of the CCFM equation (0) (continuous curve) and of the approximate 
form of this equation corresponding to setting Q{Q — q) in place of Q{Q — qx/x') 
and Aji = 1 (dotted curve). We also show the values of F2 obtained from coUinear 



factorization using the MRS (A') [^ and GRV [^ partons (the figure is taken from 
ref. m). 



3. The ladder diagram with quark (antiquark) exchange along the chain. 
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4. The effective slope X{x,Q^) defined by the formula (PB]) of the solutions of the 
equations (|^) (dashed curve) and (|35|) (solid curve). The dashed-dotted curve 
corresponds to the approximation of eq. (0) in which only the first integral on 
right hand side is retained. The dotted curve represents the slope of the solution of 
eq. ( P^D in which the inhomogeneous term is set proportional to The slopes 

are plotted as functions of x for fixed = lOGeV^. 
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Fig. 2 
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